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Abstract 



This paper presents a novel analysis of the action/energy density distribution around 
a static quark- antiquark pair in SU(3) lattice quantum chromodynamics. Using 
a , the Fourier transformation of the link variable, we remove the high-momentum 

_ gluon and extract the flux-tube component from the action/energy density. When 

Q ' the high- momentum gluon is removed, the statistical fluctuation is drastically sup- 

. pressed, and the singularities from the quark self-energy disappear. The obtained 

flux-tube component is broadly distributed around the line connecting the quark 
CN ■ and the antiquark. 
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1 Introduction 



Quark confinement is one of tlie most significant phenomena in quantum chro- 
modynamics (QCD). Due to the non-perturbative and non-Abelian nature of 
QCD, analytical derivation of quark confinement from the QCD Lagrangian 
has not yet been achieved. While it has long been known that quark confine- 
ment is well described by a "string" or "flux tube" , its mechanism remains an 
unsolved problem [1]. 

Lattice QCD provides us with beneficial knowledge about quark confinement. 
In lattice QCD, the formation of the flux tube is visualized by the analy- 
sis of the action/energy density distribution around a static quark- antiquark 
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pair [2|3Pf5|6f7f8] ■ At the positions of the quark and the antiquark, the ac- 
tion/energy density is strongly enhanced due to the divergence of the quark 
self-energy. In addition to such large perturbative contributions, a string-like 
structure is formed between the quark and the antiquark. This is direct evi- 
dence of the flux-tube formation. 

In recent works of lattice QCD, the static quark-antiquark potential is ana- 
lyzed by the Fourier transformation of the link variable [9|T0] . These works 
clarify that quark confinement originates from an infrared gluon below 1.5 
GeV in the Landau gauge. When this infrared gluon is removed, the confine- 
ment potential completely disappears. Conversely, when a high-momentum 
gluon above 1.5 GeV is removed, the short-range Coulomb potential disap- 
pears, and the quark-antiquark potential becomes a purely linear confinement 
potential. Thus, by restricting the gluon field to the low-momentum region be- 
low 1.5 GeV, we can extract the essential contribution to color confinement. 

In this letter, we apply this type of analysis by the Fourier transformation to 
the calculation of the action/energy density distribution. Because the flux-tube 
contribution to the action/energy density is rather smaller than the pertur- 
bative contribution, it is difficult to observe the flux-tube structure in tradi- 
tional approaches. For example, the endpoints of the flux tube are completely 
hidden by the perturbative singularities. If the high-momentum gluon is re- 
moved, the unnecessary perturbative contribution would disappear while the 
flux-tube contribution would remain unchanged. Therefore, by removing the 
high-momentum gluon, we can clearly observe the flux tube. In other words, 
we can extract only the flux-tube component from the action/energy density 
distribution. 



2 Action density and energy density 

In lattice QCD, the action density is deflned as 

p{x) = /3 E |l - ^ReTrf/^,(x)| , (1) 

where U^^{x) is the plaquette variable and (3 = 2Nc/g'^. The spatial distribu- 
tion of the action density around a static quark-antiquark pair is obtained by 
measuring p{x) around the Wilson loop at a certain time slice. Its expectation 
value is given by 

, , (p(x)W(R,T)) , , 

{pix))w ^ - (2) 
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Fig. 1. The Wilson loop W{R,T) and the three-dimensional plane where the action 
density {p{x))w is measured. The origin of the coordinate is placed in the center of 
the Wilson loop. 

where W{R, T) is the value of the Wilson loop with the size RxT. We consider 
only the case that the spatial size R and the temporal size T are even numbers 
in the lattice unit. The schematic figure is illustrated in Fig. [TJ The origin of 
the four-dimensional coordinate is placed in the center of the Wilson loop. 

The action density distribution {p{x))w is independent of the time slice if the 
ground-state component is suitably dominated. The temporal size T should 
be sufficiently large to extract the ground-state distribution. In the following, 
{p{x))w is measured at the central time slice of the Wilson loop {t = 0). 

The energy density e{x) is calculated identically to the action density by 
changing the relative sign between the spatial plaquettes and the temporal 
plaquettes. Because these distributions around a static quark- ant iquark pair 
are qualitatively similar except for the overall values, we mainly discuss the 
action density in the following. 



3 Momentum cutoff for the link variable 

By removing the high-momentum gluon, we investigated the behavior of the 
low-momentum gluon, which contributes to the confinement potential and the 
fiux tube. The lattice framework to remove momentum components of the 
gluon field is introduced in recent works piTO] . We used the three-dimensional 
version of this framework, i.e., we treated the spatial three-momentum of the 
gluon field. Here, we only briefiy explain the procedure. The procedure is 
formulated as the following five steps. 

Step 1. The SU(3) link variable U^{x) is generated by Monte Carlo simulation. 
Here, the link variable must be fixed with a certain gauge. In this letter, we 
choose the Landau gauge, which has a direct connection between the link 
variable and the gauge field. 
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Step 2. The momentum-space link variable Uf^{t,p) is given as 



Un{t,p) = J^Yl ^m(^) exp(ip • x), (3) 

X 

where L is the number of lattice sites in the spatial direction. The momentum 
space is also an lattice, and its lattice spacing is given by ttp = 27i/La. 

Step 3. The high- momentum component of U^j,{t,p) is removed above the ul- 
traviolet cutoff Auv- The momentum-space link variable with the ultraviolet 
cutoff is defined as 

i UJt,p) dpi < Auv) 
U^it,P)-{ (4) 
I (|p1>Auv). 



Step 4. The inverse Fourier transformation is performed as 

W = E^^(^'P)expHp-f). (5) 

Note that Ul^{x) is not an SU(3) matrix in general. To obtain the SU(3) link 
variable, U'^{x) is projected onto an SU(3) element U^{x). Such a projection 
is realized by maximizing the quantity 

ReT.[{Ui{x)}%{x)]. (6) 



Step 5. The expectation value of the operator O is calculated from the link 
variable U^{x) instead of C/^(x); i.e., {0[U^]) instead of {0[U]). 

Using this framework, we analyzed the low-momentum part of the action/energy 
density distribution. This is one type of Fourier analysis of a spatial distri- 
bution. However, unlike ordinary Fourier analysis, its Fourier component is 
not that of the spatial distribution itself but rather of the gluon field. In this 
sense, our analysis is a physical extraction rather than a numerical technique. 

We comment on the gauge fixing in Step 1. In general, because the gauge 
transformation is nonlocal in momentum space, the momentum region of the 
gauge field is a gauge- dependent concept. Therefore, the gauge fixing is nec- 
essary to remove a part of the momentum region. In this letter, we show the 
numerical results of the Landau gauge. It is noteworthy, however, that one 
can easily analyze another gauge in the same way. 
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It is worth noting that, as a by-product, our approach strongly suppresses 
the statistical noise of the gluon distribution. Due to the large fluctuation 
of the ultraviolet gluon, precise observation of the flux tube is not easy in 
lattice QCD. As shown in Sec. 5.3, the statistical noise is indeed suppressed 
by removing the ultraviolet gluon, and the gluon distribution becomes clear. 



4 Simulation setup 

We performed quenched lattice QCD simulation. The gauge action is the 
SU(3) plaquette action with (3 = 6.0. The corresponding lattice spacing a 
is about 0.10 fm, which is set so as to reproduce the string tension a = 0.89 
GeV/fm. The lattice volume is mainly 16"^. The configuration number is 500, 
and the statistical error is estimated by the jackknife method. For statistical 
improvement, averages over the translationally invariant quantities are taken, 
if possible. For example, {p{x)W{R,T)) is given by the convolution sum as 

{pix)WiR, T)) = {^Y: p{x + s)W{R, T, s)), (7) 

where s is the position of the parallel-translated Wilson loop. 

To enhance the ground-state component, the APE smearing method is applied 
to the spatial link variables of the Wilson loop [11]. In the calculation of the 
gluon distribution, it is necessary that the ground-state component is suffi- 
ciently dominant in the time slice where the gluon distribution is measured. 
Due to the smearing method, the Wilson loop is rapidly dominated by the 
ground-state component. In fact, the gluon distribution is found to be almost 
independent of T in the range of T > 4a in the present calculation. 



5 Numerical results 

5.1 Quark- antiquark potential 

First, using the framework explained in Sec. 3, we demonstrate how the quark- 
antiquark potential is changed by the ultraviolet cutoff. The original quark- 
antiquark potential is known to be expressed as 

V{R) = aR-^ + C, (8) 
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Fig. 2. The quark-antiquark potential V{R) with the ultraviolet cutoff Auv- The 
"No Cut" is the original lattice QCD result. 

where R is the interquark distance. The confinement potential is an infrared 
property. On the other hand, the one-gluon-exchange Coulomb potential and 
the constant term, which originates mainly from the regularization for the 
short-range singularity, are ultraviolet properties. 

In Fig. [21 the quark-antiquark potential with the ultraviolet cutoff Auv is 
shown. The result of the original lattice QCD is also shown in the figure. 
When the high-momentum gluon is removed, the Coulomb potential and the 
constant term disappear because these terms are ultraviolet properties. On the 
other hand, the confinement potential is almost unchanged up to Auv = 1-5 
GeV. At Auv = 1-0 GeV, the slope of the confinement potential suddenly 
decreases. In Table [H we list the asymptotic string tension (Tasym, which is 
estimated by fitting the quark-antiquark potential in i? > 0.4 fm with a linear 
function (Xasym-R + coust. In Auv > 1-5 GeV, the asymptotic string tension is 
almost equal to the original value. 

This behavior of the quark-antiquark potential is consistent with the case of 
the four-dimensional cutoff piTOj . In the Landau gauge, the low- momentum 
gluon below 1.5 GeV is relevant for quark confinement from both the three- 
dimensional and the four-dimensional viewpoints. 



5.2 Action density in vacuum 



Before the action density distribution with a static quark-antiquark pair, we 
calculated the action density without color sources, i.e., in vacuum. The vac- 
uum action density relates to the gluon condensate in continuum QCD, which 
produces the trace anomaly [T2|[T3j. In the naive continuum limit, the action 
density corresponds to the gluon condensate; however, at finite lattice spacing. 
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Fig. 3. The action density distribution {p{x))w around a static quark-antiquark 
pair. The separation between the quark and antiquark is -R = 0.6 fm. "x5"s mean 
that the lower figures are five times enlarged in the vertical direction compared to 
the upper figures. The "No Cut" (the upper left) is the original lattice QCD result 
without the ultraviolet cutoff, and its statistical error is relatively large. 

it is dominated by perturbative corrections rather than the non-perturbative 
gluon condensate. In short, a large part of the vacuum action density on the 
lattice is perturbatively generated. 

The vacuum action density {p{x)) with the ultraviolet cutoff is hsted in Table 
[H In the table, the statistical error is omitted because it is negligibly small. 
When the high-momentum gluon is removed, the vacuum action density dras- 
tically decreases. For example, at Auv = 1-5 GeV, the vacuum action density 
is reduced to about 1 % compared to the original value. The remaining small 
component would lead to non-perturbative properties of QCD vacuum. 

Table 1 

The ultraviolet cutoff Auv; the asymptotic string tension o"asym of the quark- 
antiquark potential, and the vacuum action density {p{x)). The statistical error 
of {p{x)) is omitted because it is negligibly small. 

Auv [GeV] dasym [GcV/fm] {p{x)) Ja=^ 



No Cut 0.89 14.51 

3.8 0.824(31) 2.57 

2.2 0.801(67) 0.60 

1.5 0.799(18) 0.20 

1.0 0.208(4) 9.5x10'^ 
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5.3 Action density with a quark- antiquark pair 

The action density distribution around a static quark-antiquark pair is mea- 
sured by Eq. (|2]). The ultraviolet cutoff is introduced to the action density 

as 

The arguments, such abbreviated for simplicity. The physical in- 

terpretation is the spatial distribution of the low-momentum gluon around a 
physical quark-antiquark pair. As another choice, one can introduce the ultra- 
violet cutoff not only for p[U] but also for iy[?7], and its result is expected to 
be qualitatively similar to that of Eq. 

In Fig. 131 we display the action density distribution {p{x))w with the ultra- 
violet cutoff Auv = 1.5 GeV, 2.2 GeV, and 3.8 GeV. The interquark distance 
between the quark and the antiquark is i? = 0.6 fm. The overall sign of {p{x))w 
is flipped in the figure, which is only a matter of definition. Because the abso- 
lute value of {p{x))w at Auv = 1-5 GeV and 2.2 GeV is small, these data in 
Fig. |3] are enlarged by a factor of five compared to the other ones. The action 
density distribution in original lattice QCD without the cutoff is also displayed 
("No Cut" in the upper left); however, its statistical error is relatively large. 
In original lattice QCD, the action density is strongly enhanced at the posi- 
tions of the quark and the antiquark. In contrast, the flux-tube structure is 
difficult to observe due to such singular peaks and the large statistical fluc- 
tuation. Both the singular peaks and the large fluctuation originate from the 
perturbative property of the action density. When the high-momentum gluon 
above 3.8 GeV is removed (the upper right), these perturbative contributions 
are drastically suppressed, and the flux-tube structure connecting the quark 
and the antiquark becomes clear. At Auv = 1-5 GeV and 2.2 GeV (the lower 
right and the lower left, respectively), the two peaks seem to disappear, and 
the action density is distributed around the origin. 

Apart from the vacuum contribution, which is translationally invariant, the 
action density distribution at Auv = 1-5 GeV is broadly distributed around the 
midpoint between the quark and the antiquark. In the calculation of the quark- 
antiquark potential, this low-momentum gluon leads to the linear confinement 
potential over the entire range of R. Therefore, it is considered that this action 
density distribution corresponds to the confining flux tube. 

In Fig. m we plot the action density distribution {p{x))w along the x-axis, 
i.e., in the longitudinal direction of the quark-antiquark separation. In origi- 
nal lattice QCD without the ultraviolet cutoff ("No Cut"), the action density 
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Fig. 4. The action density distribution {p(x))w along the x-axis. The interquark 
distance R is 0.6 fm. The quark and the antiquark are located on x = 0.3 fm and 
X = —0.3 fm, respectively. The "No Cut" is the original lattice QCD result. 
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Fig. 5. The action density distribution {p{x))w and the energy density distribution 
{e{x))w along the x-axis. The interquark distance R is 0.6 fm. The quark and the 
antiquark are located on x = 0.3 fm and x = —0.3 fm, respectively. 



has two self-energy peaks at the positions of the quark and the antiquark, 
X = 0.3 fm and x = —0.3 fm, respectively. When the ultraviolet cutoff is 
introduced, these self-energy peaks are drastically suppressed. Moreover, the 
absolute values of the action density and the statistical fluctuation become 
small, as in the case of the vacuum action density. The results at Auv = 1-5 
GeV and 2.2 GeV are also shown in Fig. [51 The action density distribution 
has a maximum at the origin, and the self-energy peaks seem to disappear. 
Although the self-energy peaks would also include a nonperturbative contri- 
bution, it is too small to distinguish from the flux tube. The endpoints of the 
flux tube are not sharp, spreading outside the positions of the quark and the 
antiquark. In the case of i? = 0.6 fm, the longitudinal shape of the flux tube 
resembles a broad mountain rather than a plateau. 
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Fig. 6. The action density distribution {p{x))w and the energy density distribution 
{e{x))w along the y-axis. The interquark distance R is 0.6 fm. The sohd hues are 
the results of fitting to a the Gaussian function po exp(— 



We also show the energy density distribution {e{x))w along the x-axis in Fig.|5l 
The absolute value of the energy density is smaller than that of the action 
density due to the cancellation between the chromoelectric contribution and 
the chromomagnetic contribution. Apart from the absolute value, the overall 
shape of the energy density distribution is similar to that of the action density 
distribution. 



To estimate the width of the flux tube, we fit the action density distribution 
along the y-axis to a the Gaussian form, po exp(— It is seen that the 
Gaussian form can well reproduce the lattice data, as shown in Fig. El The 
best-fit width parameter 5 is 0.31±0.01 fm at Auv = 2.2 GeV, and 0.35±0.01 
fm at Auv = 1-5 GeV. These values are comparable to the flux-tube width of 
earlier works in the standard lattice QCD [if5|6] . 

Next, we analyze how the flux-tube shape depends on the interquark distance 
R. As shown in Fig. |3l the longitudinal length and the transverse width are 
almost the same at i? = 0.6 fm and Auv = 1-5 GeV. The overall shape seems 
to be isotropic and far from a "string" or "tube" . This is because the transverse 
width of the flux tube is fairly large. To approach the gluon distribution to 
the tube-like shape, the interquark distance must be enlarged |SE]- Here, we 
use a 32(x-axis)xl6^ lattice instead of a 16^ lattice and extend the interquark 
distance to R = 1.0 fm. As shown in Fig. [TJ the distribution is stretched in 
the longitudinal direction, and the longitudinal length becomes larger than 
the transverse width. Compared to the result oi R = 0.6 fm, the flux tube at 
i? = 1.0 fm approaches a broad tube- like shape. 
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Fig. 7. The action density distribution {p{x))w for the interquark distance i? = 1.0 
fm. This calculation is performed on a 32x16^ lattice. 

6 Summary 

Using the Fourier transformation of the link variable in SU(3) lattice QCD, 
we have extracted the fiux-tube component from the action/energy density 
distribution. By removing the high-momentum gluon above 1.5 GeV in the 
Landau gauge, we can eliminate the unnecessary perturbative contribution, 
such as the singular peaks at the positions of the quark and the antiquark. As a 
by-product, the ultraviolet statistical fluctuation is also drastically suppressed. 

The resultant flux-tube component is broadly, almost isotropically, distributed 
when the interquark distance is not large, as shown in Fig. [31 When the in- 
terquark distance becomes larger, the fiux-tube component is stretched in the 
longitudinal direction and approaches a broad tube-like shape, as shown in 
Fig. [71 These distributions are the essential shapes of the confining flux tube. 

Finally, we comment on the Liischer term, which is a — 7r/(12i?) correction to 
the quark-antiquark potential due to the string fluctuation p3fT5] . Its func- 
tional form is similar to the perturbative Coulomb potential, but its origin 
is nonperturbative. While the short-range Coulomb potential vanishes due to 
the ultraviolet cutoff, as shown in Sec. 5.1, it is nontrivial whether or not the 
long-range Liischer term remains. It is interesting to analyze the Liischer term 
in this framework. 
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